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We develop a non-equilibrium quantum model of Bose-Einstein condensation and lasing of photons 
in a dye filled microcavity. We examine in detail the nature of the thermalization process induced 
by absorption and emission of photons by the dye molecules, and investigate when the photons 
are able to reach a thermal equilibrium Bose-Einstein distribution. At low temperatures, or large 
cavity losses, the absorption and emission rates are too small to allow the photons to reach thermal 
equilibrium and the behavior becomes more like that of a conventional laser. 



Bose-Einstein condensation (BEC) has been observed 
in a wide variety of systems, from ultra-cold atomic 
gases [H, m to quasi-particles in solid state systems such 
as polaritons [3|-|6|, excitons [7] and magnons [Sj. Re- 
cently experiments have also shown convincing evidence 
of a Bose-Einstein distribution, and macroscopic occu- 
pation for a gas of photons confined in a dye-filled op- 
tical microcavity |9Ml2]. In these experiments, the ther- 
mal equilibrium distribution of photons arises because 
of phonon-dressing of the absorption and emission by 
the dye molecules, and the rapid thermalisation of vi- 
brational modes of the dye molecules by their collisions 
with the solvent. This leads to the accumulation of low 
energy photons, closely following a Bose-Einstein distri- 
bution, as is clearly seen experimentally [iq| . 

Such a system is very closely related to a dye laser [l3| , 
but differs in the near-thermal emission that is observed 
below and near threshold. In the context of polariton 
condensation 0-0 there has been much debate [alii 
about the extent to which the lack of thermal equilib- 
rium in polariton experiments means the system should 
be called a condensate or a laser. However calculations 
for polaritons using various methods, from quantum ki- 
netics H, [ill to Schwinger-Keldysh path integrals [l8| . 
have shown that there is a relatively smooth crossover 
between behavior typical of a laser, and that typical of 
an equilibrium condensate. 

Both lasers and condensates involve a spontaneous 
phase symmetry breaking, and a transition to a macro- 
scopically occupied quantum state, and so, as has long 
been recognized [I^, they share many features. However, 
recent experiments on polariton and photon condensates 
also show what important differences arise, for example, 
showing how thermalisation may lead to lasing without 
inversion. In the context of polariton condensates, this is 
closely connected with the possibility of lasing in strong 
coupling [l8|. The coupling in the photon condensate, 
however, remains weak, and the thermal Bose-Einstein 
distribution seen in this system can be understood as 
arising from the combination of asymmetry between ab- 
sorption and emission and the re-trapping of fluorescence. 
A theoretical description of how asymmetric absorption 
and emission could be engineered in a circuit QED setting 
was recently discussed by Marthaler et al [20[. In this 



Letter, we extend such an approach to understand lasing 
without inversion, thermalisation, and the limits in which 
equilibrium behavior is seen for the dye-filled optical cav- 
ity. We find that over a range of parameters (including 
those used in the experiments) the system shows BEC- 
like behavior. However, significant deviations from this 
occur if the cavity losses are made larger, and a crossover 
toward more standard lasing behavior is observed. 

Some theoretical work has attempted to produce mod- 
els of this system from the point of view of equilibrium 
statistical mechanics [2l|, [22| , while other work has exam- 
ined the emergence of phase coherence in a BEC where 
particles interact through an intermediate medium (23| . 
We aim instead to provide a general non- equilibrium 
framework for understanding the steady state properties 
of the photons, taking into account the pump and decay 
processes. This allows one to understand how these com- 
pete with thermalization process, and control when the 
system behaves like a laser or like a condensate. 

A schematic diagram of our model is shown in Fig. [TJ 
This consists of a photon modes, with creation operators 
a]^, confined in the optical cavity coupled to a single elec- 
tronic transition of the dye molecules. Dye molecules are 
represented as two- level systems, written as Pauli matri- 
ces (Ji with splitting A between ground and excited lev- 
els. These levels are dressed by ladders of rovibrational 
states, which can be thought of as an on-site phonon [l3| . 
described by operators bi^b]. The level scheme is shown 
in the insert to Fig. [H The Hamiltonian is thus: 

m i 

^dYl (^^"^^ + ^L^7) ' (1) 

As in the experiment [lo|, we consider photon modes in 
a two-dimensional harmonic potential (arising from the 
curvature of the mirrors). We therefore take regularly 
spaced oscillator levels ujm = ojq -\- me, having a degen- 
eracy gm given by ^^n = ^ + 1- The lowest frequency 
uo is the "cavity cutoff" . In the following we quote fre- 
quencies relative to the molecular splitting A and thus 
introduce Sm = — A. If in equilibrium, Bose con- 
densation would lead to a macroscopic occupation of the 
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FIG. 1. (Color online) Cartoon of the system showing the 
decay processes included in Eq. ([2]). The zoomed in view 
shows the energy level structure of the dye molecules. The 
graph shows the characteristic behavior of the emission rate, 
r{—6m) (dashed line) and the absorption rate, r{Sm) (solid 
line) described in the main text. The vertical (red) lines show 
the typical spacing of the photon modes confined in the cavity. 



lowest mode. As discussed later, the extent of thermal- 
isation is affected by the detuning of the cavity cutoff, 
Sq. Since the light-matter coupling is small compared 
to optical frequencies, we assume a Jaynes-Cummings 
coupling, with frequency independent coupling strength 
g. The vibrational mode spacing is and the interac- 
tion between electronic and vibrational states is given by 
the dimensionless Huang- Rhys parameter which char- 
acterizes the difference in phonon displacement between 
the ground and excited states. The parameter values we 
use, corresponding to the experiment [10], are given in 
figure captions. 

To model the open system, we must add additional loss 
processes and external pumping. We include the loss of 
cavity photons with rate hz^ assumed independent of the 
photon frequency, and a rate T^ describing fluorescence 
of the dye molecules into non-cavity modes. To balance 
these losses we include a pumping term with rate T^. 
These processes may all be described by standard Marko- 
vian Lindblad terms, as there is no significant thermal 
occupation of relevant photon modes outside the cavity. 
The localized vibrational modes also undergo incoherent 
relaxation, due to scattering off of the solvent molecules. 
This is modeled as a relaxation rate 7 toward a ther- 
mal equilibrium state at temperature T. These processes 
cannot be described by Markovian loss rates, as this can- 
not describe thermalisation of the radiation 24, 25|. Be- 
low we describe an alternate approach to include these 
processes. For a good cavity, where equilibrium may be 
reached, 1/k is the longest timescale. 



If the coupHng to phonons, is reasonably strong, 
then multiphonon effects will be important in describ- 
ing the thermalisation processes. These can be cap- 
tured by making a polaron transformation H WHU 
where U = ex.p[^.^/Sa^{bi — 6^)]. Since the coupling 
of molecules to the optical modes is weak, we then treat 
the dynamics perturbatively in g while keeping all orders 
of S. Working in the interaction picture, and expanding 
the Liouville equation to second order in ^, one may then 
trace out the degrees of freedom associated with the vi- 
brational mode and its damping. The resulting equation 
of motion then contains Lindblad terms which cause si- 
multaneous transitions in both the photon field and the 
dressed molecule [20|. These processes then describe the 
emission (absorption) of photons into (from) the cavity, 
as shown schematically in Fig. [TJ Including all processes, 
the resulting master equation for the photon-molecule 
system is: 

p = -i[Ho,p] - ^ If £[a™] + ^C[at] + ^C[a-] 

1,171 ^ 

+ 2 ^K^i ] + I p. (2) 

Here the system Hamiltonian is Hq = ^rndtidm and 
C[X]p= {XtX, p} - 2XpX^ is the usual Lindblad term. 
The phonon assisted emission and absorption rates are 
given by the function [2o|, 

/CO 
dt/(t)e-(r^+r^^)l*l/2e-'"*, (3) 
-00 

evaluated at the detuning of the relevant photon mode. 
Here f{t) is a correlation function of polaron opera- 
tors [ifl, S, [13 given by. 



f{t) = exp 



2^7 
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dv — 
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(4) 

An illustration of these rates, as a function of detuning, 
is shown in Fig. [H We note that for frequencies where 
the rates r^,r| can be ignored in Eq. (|3]), the vibration 
induced emission and absorption rates are related by a 
Boltzmann factor V(S) = e-^^r(-^)[2H with /3 corre- 
sponding to the phonon (solvent) temperature, thus sat- 
isfying the Kennard-Stepanov relation between absorp- 
tion and emission rates [^]. At large frequencies r(a;) 
ceases to obey this relation because the incoherent pump- 
ing process corresponds to coupling to a white noise (i.e. 
infinite temperature) bath [25|. 

Marthaler et al. [lO] considered this kind of Master 
equation as a route to lasing without inversion in circuit 
QED. The same mechanism they proposed also applies 
for the photon condensate, allowing coherent emission far 
below inversion: if r{—dm) > r(^m) then the asymme- 
try in emission and absorption induced by thermalisation 
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with the phonons allows net gain without inversion. For 
lasing to occur significantly below the inversion point we 
require Sq <C — T so that the asymmetry of the absorp- 
tion and emission rates is sufficiently large. As we will 
discuss below, the same conditions lead to thermalized 
lasing, as long as the relevant emission rates T{—Sm) are 
large enough to overcome the losses from the cavity. 

We can use the master equation, Eq. (|2]), to derive a 
semi-classical rate equation for the population of each 
photon mode which, after adiabatically eliminating the 

molecular degrees of freedom, is: 

where we define = + T^rn9m^{^m)nm and f; 
r^ + Zlm ^mr(-(5m)(n^ + l). We can then use the steady 
state of this expression in combination with the rates 
from Eq. ([3]) to calculate the photon population in each 
mode, gm'^m- Note that the Um + 1 term in the emission 
process corresponds to the trapping of spontaneous fluo- 
rescence from the dye. Along with the Kennard-Stepanov 
relation discussed above, this means that in the equilib- 
rium limit tv^r^^r^ 0, but r^/F^ finite, the below 
threshold solution to this equation is simply an equilib- 
rium Bose-Einstein distribution with an effective chemi- 
cal potential given by /ieff = T ln(r^/r|) which condenses 
when /ieff = ^0- 

In Fig. [2fa) we present results for parameter values 
(given in the caption) typical of those in the experiments 
of Refs. (ol-llH. We fit Bose-Einstein distributions to the 
data by tuning the chemical potential so that the tail 
matches the numerical results. At pump powers far be- 
low threshold the occupation of the modes is, to a very 
good approximation, given by a Boltzmann distribution 
(in this limit the /ieff <^ so the distribution is classical). 
For pump strengths above the condensation threshold 
there is a macroscopic occupation of the lowest energy 
mode. The chemical potential approaches this mode, 
producing a condensate. If the pump is increased even 
further beyond this point (so that more modes are above 
threshold) the only effect is that the occupation of the 
lowest energy mode increases further, the shape of the 
tail does not change. 

The profile changes dramatically when the cavity loss 
rate is increased, as shown in Fig. [2fb). In this case 
the losses for the lowest energy modes exceeds the gain, 
controlled by the rate F(— (5o). The loss dominates these 
modes, thus preventing the low energy modes from reach- 
ing thermal equilibrium. As can be clearly seen in the 
figure, the modes with higher energy (i.e. Sm — 0) still 
match the Bose-Einstein distributions well. At suffi- 
ciently strong pumping there is a threshold (at much 
higher power than would be required if in thermal equi- 
librium) above which we find a macroscopic peak ap- 
pears in an excited mode of the cavity. The mode that 
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FIG. 2. (Color online) Mode populations gmrim vs detuning 
5m for various pump strengths. Stars are results of the non- 
equilibrium model, and lines show Bose-Einstein distributions 
fitted to the tail of the numerical results. Insets show the 
pump powers, F^/F^ = g/^Meff^ (dashed lines) compared to the 
ratio F(— ^rn)/F(^rn) — e^*^"^ (solid, black line) which gives a 
good approximation to the threshold. Panel (a) corresponds 
to experimental losses, = lOMHz. Panel (b) shows k — 
5GHz where losses prevent thermalisation. Other parameters 
are: 7 = lOOTHz, F^ = IGHz, S = 0.5, Q = ITHz, N = 10^\ 
g = O.lGHz, T = 300K, ^0 = -200THz and the mode spacing 
e = lOTHz. 



condenses is determined by the lowest mode such that 
F(— (5m) is large enough to overcome the losses. For the 
parameters of Fig. EJb) the thermalisation still plays a 
role in the rates of emission and absorption, so that it 
is not the mode with peak emission rate (near Sm = 0) 
which condenses. However, at yet higher decay rates, 
the behavior crosses over toward such "standard" laser 
behavior, and all thermal properties are lost. Similarly, 
as T ^ the emission and absorption spectra become 
a narrow Lorentzian peak at (5m = 0, and condensation 
moves to the center of the gain peak. 

Since the origin of the destruction of thermalization is 
the competition between loss and emission rate F(— (5o), 
it is also clear that lowering the cavity cutoff (making 
So more negative) has a similar effect. Equilibrium be- 
havior can only be seen when the cavity cutoff is suffi- 
ciently close to the molecular frequency, which ensures 
that F(— ^0) is sufficiently large. 

In order to explore the degree of thermalisation as a 
function of temperature and loss rates, we next consider 
the behavior at threshold. We examine two different 
aspects; the threshold pump power (a measure typical 
when considering lasing), and the total photon density 
at threshold (a measure typical for an equilibrium con- 
densation transition). Fig. [3fa) shows the total number 
of photons in the cavity, Atot = Xlm 

gm^ni') as a func- 
tion of the pump rate, F^, at various temperatures. As 
expected, increasing the temperature of the phonons re- 
duces the asymmetry between T{Sm) and F(— (5^), and 
thus both the pump power and total density at thresh- 
old increase. To explore temperature dependence we 
identify the threshold as the lowest pump power where 
max(nm) > 2 — note that this maximally occupied mode 
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FIG. 3. (Color online) Panel (a): total number of photons 
in the cavity A^tot vs pump power. The dashed vertical lines 
show the threshold (defined by max(nm) = 2). Panel (b): 
threshold pump power vs temperature, for various cavity loss 
rates as indicated. The dashed (black) lines show the same 
threshold calculated for the equilibrium theory. Panel (c): 
population A^tot at threshold vs temperature for the same loss 
rates shown in (b) . The dashed (black) line is the equilibrium 
prediction of critical density accounting for finite mode spac- 
ing; the dash-dot (black) line is the equilibrium continuum 
limit discussed in the text. All other parameters are the same 
as in Fig. [2] 

is not necessarily the lowest energy mode. 

Figure ^h) shows the threshold pump power vs tem- 
perature at various cavity loss rates. As T decreases the 
absorption and emission rates at a given detuning Sm de- 
crease. At low temperatures this again means that the 
loss rates exceed the gain for the lowest frequencies, and 
thermalisation breaks down. The temperature at which 
this happens increases with increasing tv. For the strongly 
lossy case, shown by the dot-dashed (green) curve [cor- 
responding to the used in in Fig. [2fb)] then the thresh- 
old is shifted across the whole temperature range shown. 
In the high temperature limit the emission and absorp- 
tion rates are symmetric, so the threshold pump strength 
eventually becomes that required to reach inversion. 

A more common description of the "threshold" for 
equilibrium Bose condensation is the temperature depen- 
dent critical number of particles, i.e. the total number of 



photons in the cavity at threshold. For particles confined 
in 2D one expects A^tot oc [30]. In Fig.[3l^c) we plot this 
critical number for the same values of as used in Fig. 
[3](b). Along side these results we also plot, as the dashed 
(black) curve the same quantity calculated assuming the 
mode profile is Bose-Einstein distributed with a chemical 
potential, ji = 5q — Tln(3/2), this ensures that the low- 
est energy mode has an occupation of 2. For comparison 
we also show the simple analytic result in the continuum 
limit [30|, A^tot = Tr^T^/Ge^ where e is the mode spacing. 
This differs from the discrete results due to the different 
definitions of threshold, the discrete mode spacing, and 
the effect of very high energy photon modes, which are 
neglected in the numerical calculations. 

For small cavity losses we see that for temperatures 
above ~ 200K the agreement between the discrete Bose 
distribution and the numerics is very close. Below this 
temperature, the mode which condensed is no longer the 
lowest energy mode, and condensation requires stronger 
pumping. This leads to an increase in the critical num- 
ber of photons at low temperatures. As n is increased 
we again see that the temperature above which the re- 
sults match a Bose distribution increases. There are also 
notable oscillations in the critical number at low temper- 
atures. These occur as the mode which gains a macro- 
scopic population jumps to higher and higher energy. 

In conclusion we have presented a simple non- 
equilibrium model which accurately describes the steady 
state properties of the dye filled cavity systems used to 
observe BEC of photons. We found that, for relevant 
parameters, our model accurately predicts the onset of 
BEC and the equilibrium dependence of pump power 
and critical photon number on temperature. If the losses 
from the cavity are increased, the temperature reduced, 
or the detuning increased compared to those used in 
the experiment, then a crossover occurs toward behav- 
ior more typical of a laser, and thermalisation is sup- 
pressed. These results show that, as for polariton con- 
densation [Til, a smooth crossover between typical laser 
behavior and equilibrium condensation can arise. In con- 
trast to the polariton case, the process described here 
involves a thermalized stimulated emission of radiation, 
rather than stimulated scattering between different pho- 
ton or polariton modes [23]. Future studies of the time 
dynamics of how coherence arises, and the thermal distri- 
bution emerges following the switch-on of the pump can 
help to clarify this behavior, and can be predicted using 
the model presented here. 
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